THE ELECTRON DENSITIES OF 
PSEUDORELATIVISTIC EIGENFUNCTIONS ARE 
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Abstract. We consider a pseudorelativistic model of atoms and 
molecules, where the kinetic energy of the electrons is given by 
\/ p 2 + m 2 — m. In this model the eigcnfunctions are generally 
not even bounded, however, we prove that the corresponding one- 
electron densities are smooth away from the nuclei. 



1. Introduction and results 

It was proved recently [H H] that the one-electron densities of atomic 
and molecular eigenstates are smooth away from the nuclei (actually, 
real analyticity was proved in [5]). The model studied was the non- 
relativistic Schrodinger operator with fixed nuclei. The proofs in [31 
Hj depend heavily on special properties of the non-relativistic kinetic 
energy operator —A. However, the strategy of large parts of the proof 
is very robust. In the present paper we generalise the result to the case 
of so-called pseudorelativistic molecules. 

We consider an iV-electron molecule with L fixed nuclei. The pseudo- 
relativistic Hamiltonian is (in units where h — c — 1) given by 

N L „ 



H„, l (R,z) = E {r fe )-£A} + 



. OC n t^jf I -* 1 2C 7 OC q I 

, . 1=1 1 ■> 1 l<i<j<N 1 Jl 

(1.1) 

where the kinetic energy T(pj) of the j'th electron is given by the 
operator 

T{jp) = a/ 'p 2 + m 2 — m = \J — A + m? — m , 
with m G [0, do) being the mass of the electron; a is the fine structure 
constant (in these units, a = e 2 , with e the unit charge). In ( 11. II) . 
R = (i?i, R 2 ,...,Rl) G R 3L , Re 7^ Rk for k ^ £, denote the positions of 
the L nuclei whose positive charges are given by Z = (Z\, Z 2 , . . . , Z L ). 



Date: November 27, 2008. 

© 2008 by the authors. This article may be reproduced in its entirety for non- 
commercial purposes. 



2 



S. FOURNAIS AND T. 0STERGAARD S0RENSEN 



The positions of the N electrons are denoted by x = (x±, X2, ■ ■ ■ , xjy) G 
R 3N where Xj denotes the position of the j'th electron in IR 3 ; Aj is 
the Laplacian with respect to Xj. We write V = (Vi, . . . , Vat) for the 
gradient operator in R 3N . In fll.il) we have omitted the nucleus-nucleus 
interaction, Yli<k rffzjfi ; since this is just an additive constant. 

The natural space for studying the operator H A r i (R, Z) is, in view 
of the Pauli Exclusion Principle, the antisymmetric spinor space, 
Af =1 L 2 (R 3 ;C 2 ), however, our results will not depend on spin and we 
do therefore not impose this antisymmetry condition. Instead we work 
on the space L 2 (R 3N ). 

We will assume that < Z e a < 2/vr for all £ e {1, ... , L}0 In this 
case we get from [2] Proposition 2.2] (see also [6] and |TT] for the case 
of Hydrogen) that the negative Coulomb potentials constitute a small 
form perturbation of the (total) kinetic energy (i.e., it is relatively form 
bounded with relative bound less than one). The electron-electron 
interactions being positive, and relative form bounded too, we get that 
the quadratic form 

N N L „ 

q(u, v) :=(u,J2 Tito) «) ~ (« , £ E T^Tr] v ) ( L2 ) 

3=1 ' j=l 1=1 3 £l 

is closed and semi-bounded. Here, (•, •) is the scalar product in L 2 {R 3N ). 
Hence, we can define the operator H = Hjv,l(R, Z) as the correspond- 
ing (unique) self-adjoint operator. It satisfies 

H\R 3N ) C P(H) C H 1/2 (R 3N ) , 

and 

q(u, v) = (u, Hf ) , ueD(H), ueH 1/2 (R 3N ). (1.3) 

Here, T>(H) denotes the operator domain of H; we denote its form 
domain by Q(H). All this follows from (the statements and proofs 
of) [91 Theorem X.17] and [TQl Theorem VIII. 15]. See [8] for further 
references on Hjv^R, Z). 

Suppose if; E L 2 (R 3N ) is an eigenfunction of H, i.e., there exists 
E G R such that 

H0 = Eip . 



The experimental value of the fine structure constant is a i=s 1/137. For this 
value of a, 2/ (not) » 87 



PSEUDORELATIVISTIC ELECTRON DENSITIES 



3 



We define the one-electron density p G L 1 (M 3 ) (associated to ip) by 

N 

3=1 
N 

= 2_\ / \ip(xx, . . . ,xn)\ 2 5(x — Xj) dxi • • ■ dxjsf . (1.4) 
J^x J*? N 

The main result of this paper is the following. 

Theorem 1.1. Let ip G L 2 (M. 3N ) be an eigenfunction of H. Let the 
associated density p be as defined in ( 11.41) . 
Then 

P eC 00 (R 3 \{R 1 ,...,R L }). (1.5) 

Remark 1.2. 

(i) Theorem 11.11 will follow from the more general abstract Theo- 
rem 12.21 below. 

(ii) We state Theorem II .11 for Coulomb interactions, but it holds for 
more general potentials. For instance, one can use the Yukawa 

potential e , with c > 0, in one or all of the two-particle in- 

\ x \ 

teractions. See Theorem [272] below for a more general statement 
of the result. 

(iii) Since we are only interested in regularity properties of p, we 
can study each of the (finitely many) terms in fll.4[) separately. 
We will restrict ourselves to proving the statement in ( 11.51) for 

Pi(x) : = / \ip(x, x 2 , ■ ■ ■ ,x N )\ 2 dx 2 ■ ■ -dx N , (1.6) 



the proof for the other terms being analogous. Furthermore, to 
simplify the presentation, we limit ourselves to the atomic case 
(L = 1, Ri = 0, Zx = Z, < Za < 2/tt). 

Notation. We denote by B°°(U) the smooth functions with bounded 
derivatives on the open set U, i.e., 

B°°(U) = {u G C°°(U) | d a u G L°°(U) for all a } . 

2. The abstract theorem 

Our main interest in this paper is the regularity of one-electron den- 
sities of pseudorelativistic atoms and molecules with Coulomb interac- 
tions, as stated in Theorem 11.11 However, our result holds in a more 
general case, which we will state here. 
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It is known that, in the case of relativistic atoms, the potential en- 
ergy is not a small operator perturbation of the kinetic energy, if the 
values of a,N, and Z become too large. (This is also the case in 
other relativistic models than the one studied here.) In this case, as 
discussed in the introduction, the Hamiltonian is only defined as the 
(unique) self-adjoint operator associated to a semi-bounded closed qua- 
dratic form. On the other hand, the pseudorelativistic kinetic energy 
has an extra, important property: It is the generator of a positivity 
preserving semigroup. 

Our abstract conditions below are thus based on the kinetic energy T 
below being the generator of a positivity preserving semigroup. This 
fact follows from the explicit formula for the integral kernel of the 
semigroup generated by T(p); see e.g. [7, 7.11(11)]. 

The Hamiltonians considered will be of the form 

H = T + V, (2.1) 
where (with p = (p x , . . . ,p N ) e IR 37V ) 

N N 

T = T(p) = T (Pj) = E y/~ A i + m2 ~ m . ( 2 - 2 ) 

3=1 3=1 
N 

V = V(x) = £V i (x i )+ £ ir,,(.r ; - .''/,). (2.3) 

3=1 l<j<k<N 

The following are the assumptions on the potential V. 
Assumption 2.1. (i) - For all j e {1, . . . , N}, 

Vj e C°°(R 3 \ {0}) n £°°(M 3 \ B{0, 1)) . 

— For all Q C {1, . . . , N}, the quadratic form on ®,- e QL 2 (IR 3 ) 
given by the multiplication operator Vq := J2jeQ Vj( x j) * s 
a small form perturbation of Tq := J2jeQ \Pj \ ■ 

(ii) For all j, k e {1, . . . , N} with j ^ k, 

— Wj t k > pointwise and Wj^{x) = Wkj(-x). 

— Wjjk G C°°{W \ {0}) n £°°(M 3 \ 5(0,'l)) . 

— Multiplication by Wj t k defines a bounded operator from 
/f 1 (R 3 ) to L 2 (IR 3 ) (by interpolation boundedness from 
F X / 2 (R 3 ) to H-^im 3 ) therefore follows). 

Under the above assumptions it is clear that H = T + V is well 
defined as the (unique) self-adjoint operator of the the corresponding 
closed and semi-bounded quadratic form (see the introduction for de- 
tails). 

The main abstract result of this paper is the following. 
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Theorem 2.2. Let m > and let T be the (total) pseudorelativistic 
kinetic energy operator 

N 

T = V~ A i +m 2 -m. (2.4) 

Lei functions 

V r .R 3 ^R, ./.:. {1 Vj, 

W^ fc :M 3 ^M, j,A;G{l,...,iV},j^, 
6e given stzcn tnat Assumption \2.1\ is satisfied, and let 

N 

V (x) = V 3 ( x j ) + Yl W i> k ( x i ~ Xk ) ' 
i=i i<j<k<N 

Let H = T+V 6e £/ie self-adjoint operator associated to the correspond- 
ing quadratic form (closed on Lf 1 / 2 (R 3Ar ) ). Let finally ip £ L 2 (IR 37V ) fre 
an eigenfunction of H and Zet p be the associated density as defined in 

Then 

peC°°(M 3 \{0}). (2.5) 

Remark 2.3. As pointed out in Remark 11.21 (i), Theorem 11.11 follows 
from Theorem 12.21 

Proof of Theorem \2.2\ 

The smoothness of p is a direct consequence of Propostion 13.11 below. 
The argument is exactly the same as the one given in [JJ Section 3] in 
the proof of [U Theorem 1.1]. We therefore omit the details. □ 

All that remains is to (state and) prove Proposition 13.11 below. 

3. The parallel differentiation 

The fact that one is allowed to differentiate the eigenfunction ip par- 
allel to the singularities of the (total) potential V is the key ingredient 
in proving the smoothness of the density p. This approach was car- 
ried out for the non-relativistic Schrodinger operator — that is, with 
T(jPj) = — Aj in (II. ip — in [31 Proposition 1] (see also [I]). We sketch 
the main ideas before giving the exact statement of the result (Propo- 
sition 13.11 below) and its proof. 

Let u e L 2 (R d ) and V G L°°(lR d ), and assume that 



Aw = Vu . 



(3.1) 
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Then ( 13. ip implies that u G H 2 (R d ) } in particular, du G L 2 (IR d ) for any 
derivative d. Assume furthermore that, for some specific directional 
derivative d a = Ylj a j^j-> a j e ^> we nave G L°°(R a! ). As just 
argued, d a u G L 2 (R d ). Then, by differentiation of (13.11) . we find that 

A(d a u) = Vd a u + (d a V)u, (3.2) 

from which it follows that, in fact, d a u G H 2 (M. d ). Moreover, the above 
argument is easily localised: If d a V G L°°(U) for some open set U C M d , 
then we can conclude that d a u G H 2 (U). 

Using this idea (and an induction argument) on the eigenvalue equa- 
tion one finds that eigenfunctions of the non-relativistic molecular Hamil- 
tonian are smooth in certain directions and on certain open sets (see 
Proposition 13.11 for a precision of the geometry, which is the same as in 
the non-relativistic case). In the molecular case the (Coulomb) poten- 
tial is not a bounded function, but one easily sees that the argument 
carries over to the case of potentials V which are a small operator per- 
turbation of the kinetic energy. 

For the pseudorelativistic operator in (II. ip this procedure does not 
work immediately, since we cannot separate the kinetic and potential 
energies: Since the potential V is only a small quadratic form pertur- 
bation of the kinetic energy T, the operator H = T + V is only given 
as a form sum. 

The idea is then to move the term Vu to the left hand side in (13. ip to 
find the following substitute for the argument above. Let the operator 
Sj be self-adjoint with operator domain (contained in) H s (M> d ), for some 
s > 1. Suppose u G L 2 (M d ) satisfies (in the weak sense) the equation 

f)u = v E L 2 (R d ) . (3.3) 

It follows that u G V(Sj) C H s (R d ). If furthermore v G F 1 (M a! ) one can 
then take a derivative in (I3.3P and use arguments as above to conclude 
that du G £>(f>) C H s (R d ). 

However, in our case it is not easy to identify the operator domain 
of H. By the definition as a form sum, we only get that H 1 (M. 3N ) C 
X?(H) C H 1 / 2 (Ei 3N ). That is, we cannot take one derivative on some- 
thing in T>(H) as explained above and still be sure to obtain a function 
in L 2 (M. 3N ). Furthermore, the relativistic kinetic energy is not local, 
so introduction of cut-off functions in the induction argument becomes 
somewhat more complicated. 

Nevertheless, the above idea of a proof and therefore the main tech- 
nical step in [4] — Proposition 13 . II below — can still be justified. That is, 
the strategy of repeatedly differentiating an equation of the form (13.31) 
in 'good' directions remains: We partially identify the operator domain 
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P(H) in order to be able to take one parallel derivative {d Xp below) on 
functions therein. 

Proposition 3.1. Let P,Q be a partition of {1, . . . , N} satisfying 

p^0, png = 0, puq = {i,...,n}. 

Define, for P, Q as above and e > 0, 
U P {e) = {On, . . . , x N ) G M 3N I \ Xj \ > e for j G P, 

\xj - x k \ > e for j G P,k G Q) . (3.4) 

Define also 

Let furthermore H be as in Theorem \2.2\ and let if) G L 2 (R 3Ar ) be an 
eigenfunction ofH, i.e., there exists EeR such that 

Hip = Eip. 

Then 

8% p if) G L\U P (e)) for all 7 G N 3 . 

Proof. Since the proof is somewhat technical we split it in a number of 
steps in order to make the structure more transparent. We first prove 
a lemma on localization. 

Lemma 3.2. Let if G B°°(R iN ) andu G P(H). Then ipu G P(H) and 

H(ipu) = ip(Hu) + Bu , (3.6) 

where B G £>(L 2 (R 3JV )) is i/je commutator [T, <£>]. 

Proo/. Notice first that (pu G Q(H) since u G Z>(H) C Q(H) = 
H l / 2 {E? N ) and multiplication by maps H S {R 3N ) into itself for all 
s G R. Let v G P(H) C Q(H), then also Tpv G Q(#), and, since 
u G £>(H) and q is symmetric (see ( II. 2p ). 

q((pu,v) = (ipu, Ht>) , q(u,Tpv) = (Hu,Tpv) . (3.7) 

Now, we can calculate on a form core (C^°(lR 3Ar )) to obtain 

q((pu,v) — q(u,Tpv) + (Bu,v) , (3.8) 

where B is the operator [T, ip] , which is bounded on L 2 (R 3Af ) since 
<p G 23°°(R 3Ar ) (see Lemma D below) . It follows from and (l3~8l 
that 

(pu, Hv) = (pHw + 5m, v) for all v G £>(H) . (3.9) 
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Since y?Hw + Bu E L 2 (R 3N ) and P(H) is dense in L 2 (R 3JV ) we deduce 
from (JS2D that y?u E P(H*) = X>(H) and that holds. This proves 
the lemma. □ 

An auxiliary operator. We introduce the following two operators: 
Hq = {T(pj) + Vj^Xj)) + Yl W j7k ( Xj - x k ) , (3.10) 

jeQ j,k€Q,j<k 

on <g> ieQ L 2 (R 3 ), and 

Hp = J2T(p 3 )+ Yl ^j,k{xj-x k ) : (3.11) 

jeP j,k€PJ<k 

on ®, eP L 2 (R 3 ). 

By Assumption 12.11 (notice that the Wj k are non-negative, and that 
T(pj) — \pj\ is a bounded operator on L 2 (R 3 )) the quadratic form defined 
by Hq is closed and bounded from below on if ^(R 3 '^'). The operator 
Hq is then defined as the (unique) self-adjoint operator associated to 
this form; see [TOl Theorem VIII. 15]. 

It follows from Lemma lA.ll in Appendix lAl that Hp is self-adjoint 
with domain if 1 (R 3 ' p '). We here used Assumption 12.11 (ii) and that 
T(p) (and therefore, J2j € pT(pj)), as mentioned earlier in this section, 
is the generator of a positivity preserving semigroup. 

Define furthermore 

H = H Q ® 1 + 1®H P 

on 

L 2 (R 37V ) ~ (® jeQ L 2 (R 3 )) <g> (® jeP L 2 (R 3 )) . 
Since Hq and Hp are bounded below, it follows from results on tensor 
products [H p. 86] that H is self-adjoint with domain 

P(H) = [P(H Q ) ®L 2 (R 3|p| )] n [L 2 (R 3|QI ) g)P(Hp)] 

C L 2 (R 3 I Q I) <g> P(Hp) = L 2 (R 3 I Q I) ® fl- x (R 3 l p l) . (3.12) 

Choose ^ G £°°(R 3 ) for j E P and V? i)fc G £°°(R 3 ) for j E P,k E Q 
(and k E P, j E Q) satisfying 

Vj = Vj on R 3 \ 5(0, e/2) and W jtk = W j>k on R 3 \ B(0, e/2) . 

This is possible by Assumption 12.11 Define finally 

H = H + /p, (3.13) 

lp(x) = ^-(ajj) + ^j'.*^ _ • 

j&P UeP,keQ)u(jeQ,keP) 
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The operator H is self-adjoint, with P(H) = T>(H), since Vj,Wj t k G 
L°°(R 3 ). We have (in the form sense) 

H = T + V (3.14) 

with 

V{x) = I P {x) + Y,V j {x j )+ ]V i-> Ar r •'•/.)• (3.15) 

jeQ (j,keP,j<k)Ll(j,keQ,j<k) 

Let q be the quadratic form associated with H. An approximation 
argument, using that C^°(R 3Ar ) is a form core for both q and q, gives 
that for u,v G H 1 / 2 ^*) with suppw C U P (e/2), 

q(u, v) = q(u, v) . (3.16) 



The parallel differentiation. Let fx, fz G C°°(R) be a partition 
of unity on R satisfying that fi is non- increasing and fx(t) = 1 for 
t < 5/4, fx(t) = for t > 2, fx + f 2 = 1. 

For e > and P C {1, . . . , N}, P / define 

¥> P , e (x):= J];/ 2 (2|x i |/e) J] / 2 (2|x J - - ar fc |/e) . (3.17) 

jeP j£P,k&Q 

Then ^ G £°°(R 3JV ) and supp^ P>e C U P (e/2). 

We will prove the following lemma, by induction in /c G NU{0}. No- 
tice that part (1) in the lemma implies that d]. p ip G L 2 (Up(e)). There- 
fore, Proposition 13.11 clearly follows once we have proved Lemma [3.31 

Lemma 3.3. For all k G PJ U {0} the following holds: 

For all e > 0, all P C {1, . . . , N} with P ^ 0, and all 7 G N 3 with 

|t| < k: 

(1) ^M6D(H)nP(H). 

(2) 1/7 = 71 + • • • + 7fc ; |7j| = 1 for all j , then 

U(d2 P (<pp,^))^Ed2 P (<pp, e iP) + d2 P [T,cp P ^ (3.18) 

k 

- E^r' +7j "{(^^)^ i+ ''' +7fe w)} • 

i=i 

Proof: We proceed by induction. 

It follows from Lemma 13.21 that the statement is correct for k = 
(in which case (13.181) reduces to (13.61) . when using that Hi/j = Eip). 

Suppose that the statement is true for some k > 0. Let 7 G N 3 with 
I7I = k, and write w 7 = d2 P (<fp, e ^>)- 
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Let ep be any of the three unit vectors in R 3Ar which define the 
directions of xp. More precisely, introduce the canonical basis for R 3Ar , 
{ej} with j G {1, • • • , N}, k G {1, 2, 3}. Then the vector ep is one of 
the three possibilities 

e P := -j= ej , A: e {1, 2, 3} . (3.19) 

V I I j £ p 

Let 9 ep = ep ■ V be the directional derivative in the direction ep, and 
define the self-adjoint operator ep • p = — iep ■ V with domain 

V(e P • p) = {/ G L 2 (R 3JV ) | 9 ep / G L 2 (R 37V )} . 

Let furthermore, for t G R, r tep be the translation operator (r tep /)(x) = 
/(x + tep). Clearly t i— > r tep defines a strongly continuous semigroup 
with generator ep • p. 

Notice that for t sufficiently small, supp r tep u 1 C Up(e/2). Since 
G ©(H) C L 2 (R 3 I Q I) (8) i? x (R 3 l p l) by the induction hypothesis, we 



know that 

S ep M 7 G L 2 (R 3N ) 

so w 7 G ©(ep • p) and 



lim -(r te pU 7 — u 7 ) = d ep u 7 , (3.20) 

in L 2 (R 3Af ). 

Let u G ©(H) and consider (Hv, d ep u^). Using (I3.20p and (I3.16p . 

we get 

(Hv, (9 e p-u 7 ) = lirni - (Ht>, r tep Uj — w 7 ) = lirnt q(u, r tep u 7 — w 7 ) 

= limt _1 q(f , r tep M 7 — w 7 ) . (3-21) 

Since the translation r tep commutes with H = Hp + Hq (see (13.101) 
and (I3.1ip ). we get that, with I P from (I3.13p . 

Hr tep = r te pH + [I P , r tep ] . 

Thus, using (13.161) 
(Hv, d ep Uj) = limt _1 q(f , r tep u 7 — m 7 ) 

= limt _1 (t>, r tep (Hu 7 ) — Hu 7 ) — (v, (d ep Ip)u 7 ) . (3.22) 



To prove that d ep u 7 G ©(H*) = ©(H) from this, it remains to show 
that when applying d ep to Hu 7 we obtain a function belonging to 
L 2 (R 3N ). Then, from (13T221 . also 

H(S ep M 7 ) = H*(S ep M 7 ) = <9 ep (Hw 7 ) - (d ep I P )u y . 



PSEUDORELATIVISTIC ELECTRON DENSITIES 11 

By (13.18p . localization and (1) from the induction hypothesis, we 
find 

H« 7 = H« 7 = Edl p ( VP ^) +c£ P [T>F, e ]V (3.23) 

k 
3=1 

We will show that when applying d ep to each term on the right side of 
( I3.23P we obtain a function belonging to L 2 (M 3N ). 

For the first term, since G P(H) C L 2 (M 3|QI ) ® #W P| ) 

by ( 13.121) and the induction hypothesis, we know that 

d ep d2 p (^) e L\R m ) . (3.24) 

For the third term, the function Ip from (13. 13[) satisfies Ip G i3°°(IR 3Ar ), 
and, as just shown, d^: p {ipp^ip) G L 2 (M. 3N ) for all \a\ < k + 1, so, by 
Leibniz' rule, 

k 

^ P (E^ + -- + ^H(^^p)^r + ^(^)}) eL 2 (R 3N ). (3.25) 

Finally, we consider the commutator term &l p [T, (pp ie ]ip in (I3.23p . 
Define ip± = (fp^/i, <f2 = 1 — <£>i- Notice that, by the definition of /i, /2, 

/ 1 (8t/e)/ 2 (2t/e) = 0. (3.26) 
By using that f\ + ji = 1 we find 

^ = E L[4-( 8 N/ e ) II 4, fc (8|^-x fc |/6), (3.27) 

where the sum is over all tuples ({s,}, {sj t k}) G {1, 2}l p l + l p l'l < 2l with at 
least one entry different from 2. Write the commutator term d2 p [T, </ } p,e]'0 
as 

^ p [T> P , e ]V = ^ p [T,^p, e ](^) +^[T^ P , e ](^) , (3.28) 
The term with (pi we write, using Leibniz' rule, as 
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SO, 



/3<7 



+ E()) [T'C^Rp^/^). 



/3<7 

By Lemma IA.2I and the induction hypothesis, we therefore see that 

d ep &l p [T,y Pie ]{y^) G L 2 (R 3N ) . (3.29) 

Now we consider the term with tp 2 in (13.281) . We will prove that also 

d ep d2 p [T, vppJ G L 2 (M. 3N ) . (3.30) 

Since T is a finite sum and y?p i<E y?2 = it suffices, up to renumbering 
of the terms, to prove that 



= d Bp 6n xp <p P Jpl + m{<fwl>) G L 2 (R 3N ) . (3.31) 



Proof of (JH3U). 
Case 1. 1 G P. 

The case P = { 1 } being immediate by Lemma IA.21 we will assume 
that Pi + 0, where P X :=P\ {1}. 

Since ypf + m commutes with multiplication operators in other 
variables, and using the support condition (13.261) . we find 



fp,eypl + mip 2 = Y[ f2(2\x j \/e) Yl fc( 2 \ x i ~ x k\/t) 
jePi jePi,k£Q 



x {/ 2 (2|xx|/e) J] f 2 (2\ Xl - x k \/e)yjpi + mf} , 

with 

/ := E f<n(*\xi\/e) J] -M 8 !^ - x k \/e) , (3.32) 

cr fceQ 

where the sum is over all a G {1, 2} 1+ '^' with a ^ (2, . . . , 2). Since at 
least one factor for each summand has to be f\ we find 

supp / C {x | min (|xi|, niin \x\ — x k \) < e/4} . 
Thus, by the triangle inequality 

supp ( J] / 2 (2N/e) J] / 2 (2|x, - a: fc |/e)/) C U Pl {t/A) . 
jePi jeP u keQ 
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Since y?p li<E /4 = 1 on Up 1 (e/4 : ) we get the identity 



( PP,e\/p1 + 'm ( P2= (<pp,e\/p 2 i + mip 2 )^pp u e/A • (3.33) 

By the induction hypothesis 

dl^ Put/ ^)eL 2 (R 3N ), (3.34) 

for all |7'| < n. Furthermore, since supp (fp ;e n supp <p 2 = 0, Lemma [A72l 
yields that 



is a bounded operator on L (R ) for all a, j3, M 
By Leibniz rule and ( 13.331) . 



E c*u**,**{9£<PP,*)\lpi + m(0> 2 )(l 

x {(l+rf)-^^/^)}, (3-35) 



01+02+03=7' 



for some constants c QljQ2ja3 . 



By definition, <9" p = J2/3<a c »,l3^i^x P f f° r some constants c^. So 
using (I3.34p and choosing d^ p = d ep d1 p and M > \*y\ + 1 in (13.351) . we 
see that 



d Bp dZ p (<pp^pl + m<p 2 ) eL 2 ^ ;x 



This finishes the proof of (13. 31 p in the case 1 G P. 
Case 2. 1 £ P. 

This case is similar but simpler than Case 1. In this case we define 
Pi = P. Arguing as previously we realize that the identity (13.331) 
remains valid. Also (I3.34j) follows from the induction hypothesis. Since 
P = Pi, we can in this case choose M = in (13.351) and get the desired 
result. This finishes the proof of ( 13.311) in the case 1 ^ P and combining 
with Case 1, we get the general result. □ 

Combining (15331) . (I3T25|) . (I3T30D . and (I3T29D . we get that 

9 ep (H« 7 )6L 2 (R 3iV ). (3.36) 

So we see from (I3T22D that for all v G X>(H), 

(Hi;, <9 ep w 7 ) = (v, <9 ep (Hw 7 ) - (d ep I P )u 7 ) . (3.37) 

From (J3T36J, (I3T3TD . and (j3HS]) we conclude that 

9 ep M 7 G £>(H*) n Z>(H*) = X>(H) n P(H) , (3.38) 
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and 

H(<9 ep u 7 ) = <9 ep (Hu 7 ) — (i9 ep /p)u 7 . (3.39) 

The equations (13.38!) an d 1 13. 39ft combine to give the statement in 
Lemma [3.31 for k + 1. 

This finishes the induction step, and by induction the statement in 
Lemma [3.31 therefore holds for all £; G N U {0}. □ 

As mentioned above, this finishes the proof of Proposition 13. 1[ □ 

Appendix A. Auxiliary results from operator theory 

In the proof of Lemma 13.11 we need the following consequence of the 
Davies-Faris Theorem (J9j Theorem X.31]). 

Lemma A.l. Suppose T > is self-adjoint with domain T>(T), and 
that T is the generator of a positivity preserving semigroup. Let V be a 
positive multiplication operator, which is bounded relative to T. Then 
H = T + V is self-adjoint on V(T). 

Proof. Choose g > such that gV is relatively bounded with respect 
to T with bound a < 1. We will prove by induction that K n = T + ngV 
is self-adjoint on T>(T) for all n G N. In order to do so, let us consider 
the following statement S(n): 

(1) K n = T + ngV is self-adjoint on V(T). 

(2) \\gV<p\\ < a\\(K n + l)<p\\ for all <p e V(T). 

(3) K n is the generator of a positivity preserving semigroup. 

Note first that S(0) is true by assumption. 

Suppose now S(n) holds true for some n > 0. By S(n) point (2), 
gV is a small operator perturbation of K n , so K n+ i = K n + gV is (by 
the Kato-Rellich Theorem [91 Theorem X.12]) self-adjoint on T>(K n ) = 
T>(T). Furthermore, using the Trotter product formula [91 Theorem 
X.51] and the induction hypothesis, it is easy to see that e~ tKn+1 is 
positivity preserving (for t > 0). Then, by the Davies-Faris Theorem 
[9j Theorem X.31], it follows that gV satisfies the bound 

\\gV<p\\ < a\\(K n+1 + 1VH for all <p G V{T) . 

Therefore S(n + 1) holds. This finishes the proof that S(n) implies 
S(n + 1) for any n > 0. 

The proof of Lemma IA.1I now follows by induction. □ 

We also state the following lemma which is used repeatedly in Sec- 
tion [3j The proof is standard and is omitted. 
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Lemma A. 2. Let \i £ £>°°(R 3Ar ) have disjoint support and let m > 
0. Then [\Jp 2 j + m, defines a bounded operator on H S (M. 3N ) for all 

set and (1 + Pj) M x[\J Pj + m > +Pj) M a bounded operator on 
L 2 {R 3N ) for allM. 

Acknowledgement. Parts of this work have been carried out at var- 
ious institutions, whose hospitality is gratefully acknowledged: Math- 
ematisches Forschungsinstitut Oberwolfach (SF, T0S), Erwin Schro- 
dinger Institute (SF, T0S), Universite Paris-Sud (T0S), and the IHES 
(T0S). SF is partially supported by a Skou Grant and a Young Elitere- 
searcher award from the Danish councils for independent research, 
a grant from the Lundbeck Foundation, and the European Research 
Council under the European Community's Seventh Framework Pro- 
gramme (FP7/2007-2013)/ERC grant agreement n° 202859. T0S is 
partially supported by The Danish Natural Science Research Council, 
under the grant 'Mathematical Physics and Partial Differential Equa- 
tions'. 

References 

[1] Anne Boutet de Monvel-Berthier and Vladimir Georgescu, Graded C* -algebras 

and many-body perturbation theory. II. The Mourre estimate, Asterisque 

(1992), no. 210, 6-7, 75-96, Methodes semi-classiques, Vol. 2 (Nantes, 1991). 
[2] Ingrid Daubechies and Elliott H. Lieb, One-Electron Relativistic Molecules 

with Coulomb Interaction, Commun. Math. Phys. 90 (1983), 497-510. 
[3] S0ren Fournais, Maria Hoffmann-Ostenhof, Thomas Hoffmann-Ostcnhof, and 

Thomas 0stergaard S0rensen, The electron density is smooth away from the 

nuclei, Comm. Math. Phys. 228 (2002), no. 3, 401-415. 
[4] , On the regularity of the density of electronic wavef unctions, Math- 
ematical results in quantum mechanics (Taxco, 2001), Contemp. Math., vol. 

307, Amer. Math. Soc, Providence, RI, 2002, pp. 143-148. 
[5] , Analyticity of the density of electronic wavef unctions, Ark. Mat. 42 

(2004), no. 1, 87-106. 
[6] Ira W. Herbst, Spectral Theory of the Operator (p 2 + m 2 ) 1 ' 12 — Ze 2 /r, Commun. 

Math. Phys. 53 (1977), 285-294. 
[7] Elliott H. Lieb and Michael Loss, Analysis, American Mathematical Society, 

Providence, RI, 1997. 
[8] Elliott H. Lieb and Horng-Tzer Yau, The Stability and Instability of Relativistic 

Matter, Commun. Math. Phys. 118 (1988), 177-213. 
[9] Michael Reed and Barry Simon, Methods of modern mathematical physics. II. 

Fourier analysis, self-adjointness, Academic Press [Harcourt Brace Jovanovich 

Publishers], New York, 1975. 
[10] , Methods of modern mathematical physics. I, second ed., Academic 

Press Inc. [Harcourt Brace Jovanovich Publishers], New York, 1980, Functional 

analysis. 



16 



S. FOURNAIS AND T. 0STERGAARD S0RENSEN 



[11] Ricardo Weder, Spectral Analysis of Pseudodifferential Operators, J. Funct. 
Anal. 20 (1975), 319-337. 

(S. Fournais) Department of Mathematical Sciences, University of 
Aarhus, Ny Munkegade, Building 1530, DK-8000 Arhus C, Denmark. 
E-mail address: fournais@imf.au.dk 

(S. Fournais on leave from) CNRS and Laboratoire de Mathematiques 
d'Orsay, Univ Paris-Sud, Orsay CEDEX, F-91405, France. 

(T. 0stergaard S0rensen) Department of Mathematical Sciences, Aal- 
borg University, Fredrik Bajers Vej 7G, DK-9220 Aalborg East, Den- 
mark. 

E-mail address: sorensen@math.aau.dk 



